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Combinatorial Ricci Flows on Surfaces 


Bennett Chow and Feng Luo 


Abstract 


We show that the analog of Hamilton’s Ricci flow in the combinatorial setting prodnces 
solntions which converge exponentially fast to Thnrston’s circle packing on surfaces. As 
a consequence, a new proof of Thurston’s existence of circle packing theorem is obtained. 
As another consequence, Ricci flow suggests a new algorithm to find circle packings. 

§1. Introduction 

1.1. For a compact surface with a Riemannian metric {X,gij), R. Hamilton in [Ha] 
introduced the 2-dimensional Ricci flow defined by the equation = —2Kgij where 
K is the Gaussian curvature of the surface. It is proved in [Ha] and [Cho] that for any 
closed surface with any initial Riemannian metric, the solution of the Ricci flow exists for 
all time, and after normalizing the solution to have a fixed area, the solution converges 
to the constant curvature metric conformal to the initial metric as time goes to infinity. 
The purpose of the paper is to study the analogous flow in the combinatorial setting. 
The simplest case of our setup is as follows. Given a triangulated surface, we assign to 
each vertex a positive radius and realize the surface by a piecewise flat cone metric with 
singularities at the vertices where the length of an edge is the sum of the two radii associated 
to the end points. We call these the (tangential type) circle-packing metrics based on the 
triangulation. We consider a natural flow (the combinatorial Ricci flow) on the space 
of all tangential type circle-packing metrics. Our results give a complete description to 
the asymptotic behavior of the solution of the Ricci flow. If we use hyperbolic geometry 
as the background metric, then for compact surfaces of negative Euler characteristic, we 
show that the combinatorial Ricci flow has solutions for all time for any initial metric and 
converges exponentially fast to the circle packing metric constructed by Thurston. 

1.2. For simplicity, we shah state our results for triangulations on compact surfaces 
without boundary. For the most general version involving cellular-decompositions, see §2, 
§5 and §6. 

Suppose A is a closed surface and T is a triangulation on X. Let V = {ui, V 2 , 
be the set of vertices in T. The set of ah edges and triangles in T are denoted by E and F. 
Unless mentioned otherwise a weight on the triangulation is a function $ : U ^ [0, 7r/2]. (In 
lemma 2.3, proposition 3.1 and proposition 4.1, we allow the weight to he ^ : E —y [0, tt)). 

Given a weighted triangulation (T, $) of a compact surface A, we consider the circle 
packing metrics on the surface A based on (T, $) dehned as follows. 

To each vertex Vi in the triangulation T, assign a positive number r^. Now realize 
each edge Cij joining Vi to Vj by a Euclidean segment of length 
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Note that, due to $(e) < 7r/2, these positive numbers Ijk, hi} satisfy triangle inequal¬ 
ities (see [Th, lemma 13.7.2]). (Also note that if the weight $ = 0, then lij = Vi + rj. This 
corresponds to the case that circle patterns come from tangency.) Thus we can realize each 
triangle AviVjVk by a Euclidean triangle of edge lengths kj, Ijk, hi- The triangle is formed 
by the centers of three circles of radii r^, rj and Vk intersecting at angles $(eij), $(ejA:) 
and ^hki) where Cnm is the edge joining Vn and Vm- This produces a flat cone metric on 
the surface X with singularities at the vertices. We call these metrics the circle-packing 
metrics (based on (T, $)). Let be the cone angle at the vertex Vi which is the sum of 
all inner angles having vertex Vi. The curvature Ki at Vi is dehned to be 27r — a^. The 
combinatorial Ricci flow is the following, 

(1.1) 5 = -A>. 

It is easy to show that for any initial circle-packing metric, the solution of (1.1) with 
the given initial value exists for all time t > 0. Our goal is to understand the asymptotic 
behavior of the solution. To this end, it is useful to consider the normalized Ricci flow, 

nv • 

( 1 . 2 ) ^^-{K,-Kav)r, 

at 

where Kav = 1/A^ the average curvature. By the Gauss-Bonnet theorem (propo¬ 

sition 3.1), Kav = 27rx(A)/A^ is time independent. The two equations (1.1) and (1.2) are 
related by a change of scales. Namely rfli) is a solution to (1.1) if and only if rfli) 

is a solution of (1.2). Furthermore, if rflt) is a solution to the normalized equation (1.2), 
the product 11^=1 ^i(t) is a constant. For this reason, when we consider the asymptotic 
behavior of the Ricci flow, it is more convenient to deal with the normalized Ricci flow 

(1.2) . Let R>o be the set of all positive real numbers. 

Definition. A solution to (1.2) is called convergent if 

(1) hm^^oo Kiit) = Kfloo) G (—cx), 27r) exists for all f, and 

(2) livat^ooTiit) = rfloo) G R>o exists for all i. 

A convergent solution is called convergent exponentially fast if there are positive con¬ 
stants Cl, C 2 so that for all time t > 0, 

\Kflt) - Kfloo) \ < cie-^^\ 


and 

Irflt) - ri(cx))| < 


Given any subset I C V, let Fi be the set of all cells in T whose vertices are in / and 
let the link of /, denoted by Lk{I), be the set of pairs (e, v) of an edge e and a vertex v so 
that (1) the end points of e are not in I and (2) the vertex u is in / and (3) e and v form 
a triangle. 
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Theorem 1.1. Suppose (T, $) is a weighted triangulation of a closed connected surface 
X. Given any initial circle-packing metric based on the weighted triangulation, the solution 
to the normalized combinatorial Ricci flow (1.2) in the Euclidean geometry with the given 
initial value exists for all time and converges if and only if for any proper subset I C V, 


(1.3) 27t\I\x{X)/\V\>- (7r-$(e))+27rx(F,). 

{e,v)eLk{I) 

Furthermore, if the solution converges, then it converges exponentially fast to the 
metric of constant curvature 2 tix{.X)/N■ 

By the work of Thurston on circle packing, condition (1.3) is equivalent to the existence 
of constant curvature circle packing metric based on the weighted triangulation on the 
surface. 

In the case that the surface has negative Euler characteristic, it is natural to replace 
Euclidean triangles by the hyperbolic ones. The combinatorial Ricci flow in the hyperbolic 
geometry is the following. Given a surface X with a weighted triangulation (T, $) where 
$(i?) C [0, 7r/2], we assign to each vertex of T a positive radius r^. We realize each triangle 
AviVjVk by a hyperbolic triangle formed by the centers of three circles of radii r^, rj and 
Tfc intersecting at angles $(ey), $(ejA:) and $(efci)- Let be the cone angle at Vi and 
Ki = 2 ti — tti be its curvature. The combinatorial Ricci how in the hyperbolic metric is, 

/. dri , 

(1.4) -—=-Kismhn. 


Theorem 1.2. Suppose (T, $) is a weighted generalized triangulation of a closed connected 
surface X of negative Euler characteristic. Given any initial choice of circle-packing met¬ 
ric, the solution to (l.f) in the hyperbolic background geometry with the given initial value 
exists for all time and converges if and only if the following two conditions hold. 

(1.5a) for any three edges ci, 62 , 63 forming a null homotopic loop in X, > tt? 

then these three edges form the boundary of a triangle of T, and 

(1.5b) for any four edges ci, 62 , 63 , 64 forming a null homotopic loop in X, ifYlt=i ^(^i) — 
2%, then these Ci’s form the boundary of the union of two adjacent triangles. 

Furthermore, if it converges, then it converges exponentially fast to a hyperbolic metric 
on X so that all vertex angles are 27 t. 

Conditions (1.5a) and (1.5b) were obtained by Thurston in his work on circle packing 
on surfaces. Indeed, according to [Th] theorem 13.7.1, these two conditions are equivalent 
to the existence of circle packing metric with curvature equal to zero at each vertex. 

The following result shows that condition (1.3) is almost the same as conditions (1.5a) 
and (1.5b) for surfaces with non-negative Euler characteristic. 
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Proposition 1.3 (a) If the Euler characteristic of the surface X is positive, then condition 
(1-3) is valid for all triangulations of X so that there are no null homotopic loops consisting 
of 3 edges ei, 62 , 63 with > tt- 

(b) If the Euler characteristic of the surface X is zero, then (1.3) is valid for all 
triangulations of X so that (1.5a) and (1.5b) hold. 

Note that if the weight is the zero function, the conditions in theorem 1.2 and propo¬ 
sition 1.3 automatically hold. 

1.3. The exponential convergence of the Ricci flow may have some applications in 
practice. Recall that in his original proof of the existence of zero curvature circle packing 
metric, Thurston gave an algorithm to hnd the circle packing. This algorithm adjusts the 
radii at vertices one at a time. It is shown in [Mo] that Thurston’s algorithm converges in 
polynomial time to the circle packing. To the best of our knowledge, the existing computer 
softwares (for instance, [St]) are based on Thurston’s algorithm. It seems Ricci flow is more 
natural and may produce a faster algorithm to hnd circle packing metrics. 

1.4. In the case that the background metric is spherical, we consider triangulations 
with zero weights. The combinatorial Ricci how is 

(1.5) ^ = -Kismri 

where we realize each triangle AviVjVk by a spherical triangle of edge lengths -|- rj, rj + 
rk, fk+ri (note that we must impose G (0, tt) and r^-l-rj-l-rA; < tt in this case). Numerical 
simulations by L.T. Cheng [Che] shows that the situation is very complicated even in the 
simplest case of triangulations of the 2-sphere, i.e., the boundary of a tetrahedron. For 
instance, if we start with the boundary of the regular spherical tetrahedron, then the 
metric will shrink to a point in hnite time. If we rescale it to have a hxed area, the cone 
angles will tend to tt (i.e., tends to the boundary of the Euclidean regular tetrahedron if 
rescaled). But in some other case, there are initial metrics so that as time goes to a hnite 
number, the metric tends to the round 2-sphere. We do not know what the asymptotic 
behavior of the solution should be in general. 

1.5. The proofs of the above theorems follow the same strategy used by Hamilton in 
his original approach in the smooth case. Namely, one hrst derives the evolution equation 
for the curvature. These equations are essentially of parabolic type. Thus the maximum 
principle can be applied to control the behaviors of the curvature. In some senses, what we 
are doing are quant ideations of Thurston’s original proof. For instance, the fact that the 
curvature evolution is a heat equation is essentially a consequence of Thurston’s lemma 
13.7.3 in [Th]. We also observed the fact that dOi/drj smh{rj) = d9j/drismh{ri) (for 
hyperbolic case) and dOi/drjrj = dOj/driri (for Euclidean case) where 6i is the inner 
angle at i-th vertex corresponding to radius circle (lemma 2.3). This was hrst observed 
in the paper of Colin de Verdiere [CV] for the zero weighted triangulation. This makes 
the Ricci how variational. Namely, it is essentially a gradient how of a convex function. 
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We observed that the degeneration condition obtained in Thnrston’s work implies that 
Palais-Smale condition holds for the convex fnnctional. Thns establishes the main resnlts. 
We remark that the convexity was hrst observed in [CV] in the case of the zero weighted 
triangnlation. Also if a triangnlation has zero weight, i.e., $ = 0, then the convergence 
part of theorem 1.2 can also be dednced from the resnlts of Thnrston [Th] and Cohn 
de Verdiere [CV]. Onr proofs are different from that of Cohn de Verdiere [CVj. Cohn de 
Verdiere nsed a resnlt in network flow to prodnce the coherent system of angles. From 
this system, he then prodnced the convex fnnction and showed that it is proper. In onr 
case, it is the maximnm principle and Palais-Smale condition which imply the exponential 
convergence. Another related work is [BSj. 

1.6. The paper is organized as follows. In §2, we collect some known resnlts concerning 
triangles. Most of the proofs will be deferred to the appendix. In §3, we derive the evolntion 
eqnation for the cnrvatnre and discnss some of its conseqnences. In §4, we establish a 
compactness theorem concerning the solntions of the Ricci flow. In §5, we prove theorem 
1.2. In §6, we prove theorem 1.1. In the appendix, we derive some basic properties of the 
triangles. 

1.7. The work is snpported in part by the NSF DMS 0203926 (at UCSD) and DMS 
0103843 (at Rntgers University). 

§2. Preliminaries on Triangles 

Let be one of the three geometries, (Enclidean), (spherical) and (hy¬ 
perbolic), i.e., a simply connected snrface with a complete Riemannian metric of cnrvatnre 
eqnal to 0, 1, and -1. For simplicity, we dehne a fnnction s{x) according the geometry 
as follows. If then s(a:) = x; if = H^, then s(a:) = sinh(a:); and if = S'^, 

then s(a:) = sin(a:). 

A fnnction is called elementary if it can be obtained from the polynomials, the expo¬ 
nential fnnctions, logarithmic fnnctions, the trigonometric fnnctions and their inverses by 
the hve operations of addition, snbtraction, mnltiplication, division and composition. 

2.1. The following lemma of Thnrston gnarantees the existence of a circle pattern in 
the geometry K^. 

Lemma 2.1. ([Th], lemma 13.7.2) For any three angles G [0,7r/2] and any three 

radii ri,rj,rk G R>o, there exists a configuration of three circles in both Euclidean and 
hyperbolic geometries, unique up to isometry, having radii ri,rj,rk and meeting in angles 
Furthermore, ifri + rj +rk < ti, then such a configuration of three circles also 
exists in the spherical geometry 3“^ and is unique up to isometry. 

Indeed, let Xi be the distance between the centers of the j-th and /c-th circles. Then 
it snfhces to verify that {xi,Xj,Xk} satishes the triangle ineqnalities that snm of two is 
larger than the third and also xi -t- Xj + Xk < ‘2 .ti in the case of 3“^. To see the hrst, we 
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have Xi < rj + r/j < Xj + Xk since < 7r/2. Furthermore, the equality Xi = Xj + Xk does 
not hold. The second inequality follows from xi + Xj x^ < 2{ri + Vj + rfc). 

Consider the triangle formed by the centers of the three circles in lemma 2.1 in K^. 
Let 9i be the inner angle at the i-th vertex corresponding to the radius circle. Then 6i 
is a smooth function of the radii r^, r^. The following is known in Thurston’s notes. 

Lemma 2.2. ([Th], lemma 13.7.3). In all three geometries or , if all angles 

^n’s are in [ 0 , 7 r/ 2 ], then we have 
(a) ddijdri < 0 , 

(h) dOijdrj > 0 fori^ j, 

(c) d{6i+6j+6k)/dri < 0 inH‘^, d{6i+6j+6k)/dri > Q in S‘^ and d{6i+6j+6k)/dri ~ 
0 in E'^. 

Furthermore, the partial derivatives dOnldvm are elementary functions in ri,rj and 
rk where n,mE {i,j,k}. 

Note part (c) follows from Thurston’s geometric argument that the area of the triangle 
increases when increases. A computational proof of le mm a 2.2 which gives the explicit 
formulas and the geometric meaning of these partial derivatives can be obtained using 
lemma Al in the appendix. See the proof of lemma 2.3. 

The next lemma was hrst obtained by Cohn de Verdiere [CV] in the case where all 
angles $n’s are zero. 

Lemma 2.3. In all three geometries E‘^, and S‘^ and for all angles ^k ^ [0, tt) 

so that the triangle of radii ri, rj and rk intersecting at these angles exists, 


d9i{drjs{rj) = d9jjdris{ri). 


We defer the proof to the appendix. In the appendix, we give two proofs. The simpler 
proof uses the mathematical software Maple VII. 

One consequence of the lemmas is the following. 

Proposition 2.4. Suppose all angles are in [ 0 , 7 r/ 2 ]. If the radii rn’s evolve in time 
according to some differential equation drn{t)/dt = —Kns{rn) where Kn is a function of 
ri, rj and rk and n = i,j, k, then the evolution of the inner angle 9i satisfies the equation, 

d9,/dt = -A,j{Kj - Ki) - A,k{Kk - Ki) - A,K,\ 

where A^s for r ^ s and A^ are positive valued elementary functions in ri, rj and rk, A is 
the curvature of the geometry K‘^. Furthermore Ars = Asr- 

Proof. By the chain rule, we have 


d9i/dt = ddijdrir'i + d9i/djr'j + d9i/drkr'k 


6 



= -ddi/dris{ri)Ki - dOi/drjs{rj)Kj - dOi/drks{rk)Kk 


= -de,ldrjs{rj){Kj - Ki) - 86,/drks{rk){Kk - K,) - D,K,. 

Here Di = ddijdrisij/) + ddi/drjsijj) + dOi/drks{rk)- By lemma 2.3, we can write Di 
as ^ % lemma 2.2 (c), the last partial 

derivative is negative in the case is positive in the case K‘^ = and is zero 

if = E"^. Dehne Ai = \Di when = 3“^ or and = 1 when K'^ = E‘^. Then 
Ai > 0. Let Amn = dOm/drnsirn) for n ^ m. Then by lemma 2.2, it is positive and is 
an elementary fnnction. Fnrthermore, the fact A^^ = A^m follows from lemma 2.3. Thns 
the resnlt follows. 

§3. The 2-Dimensional Combinatorial Ricci Flow 

3.1. Basic setup and the Gauss-Bonnet formula 

We set np the Ricci flow in the most general setting as follows. 

Snppose W is a connected compact snrface with possibly non-empty bonndary. A 
generalized triangulation T of A is a cell-division of X into cells so that the lifting of the 
cell-division in the nniversal cover is a triangnlation. In particnlar, each 2-cell (triangle) 
lifts to an embedding in the nniversal cover and there are no vertices of degree 1 or 2. Also 
there are no nnll homotopic loops formed by at most 2 edges. Fnrthermore, given three 
vertices, there is at most one triangle having them as vertices. The set of all vertices (0- 
cells), edges (1-cells) and triangles (2-cells) in T are denoted by V, E, E. Unless mentioned 
otherwise, a weight on the generalized triangnlation T is a f un ction $ : i? —[0, 7r/2]. The 
pair (T, $) is called a weighted generalized triangnlation. 

This dehnition of generalized triangnlation is slightly different from Theorem 13.7.1 
in Thnrston’s note. In fact, we were not able to rnle ont the degree 2 vertices from 
Thnrston’s dehnition. On the other hand, it is clear that if a triangnlation is realized by 
a circle packing, there cannot be any degree two vertices. 

Fix a backgronnd geometry K^. Given a compact snrface X together with a weighted 
generalized triangnlation (T, $), we assign to each vertex Vi E V = {ni, ...,nAr} a positive 
nnmber r^. If the geometry is spherical 3“^, we assnme that rj rk < holds for 
all triples {r^, r^, r^} so that the corresponding vertices form a triangle in T. 

For each triangle AviVjVk in E having edges e^j, ejk and Cki, nsing lemma 2.1 we can 
realize the triangle as the one formed by the centers of the three circles of radii ri,rj,rk 
meeting in angles ^{cij) , $(ejfc) and $(efci) in geometry K^. This prodnces a Riemannian 
metric on the snrface X of constant cnrvatnre having singnlarities at the vertices. Let 
be the snm of all inner angles having vertex Vi in triangles. Then the curvature at Vi is 
dehned to be Ki = 2 ti — ai when the vertex Vi is in the interior of X and Ki = tt — when 
Vi is in the bonndary. 

Onr hrst observation is the following. 

Proposition 3.1. (Ganss-Bonnet) Let A = —1,0,1 he the curvature of one of the three 
geometries H‘^,E‘^,3^ used in the construction. If the circle packing metric based on a 
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weighted generalized triangulation (T, $) exists where $ : i? —[0,7r), then the total curva¬ 
ture 

N 

^ K, = 27rxiX) - XArea{X). 

i=l 

Note that the three terms in the Gauss-Bonnet formula are additive with respect to 
gluing along surface boundary. Thus by doubling the surface X across the boundary if 
X has non-empty boundary, it suffices to show the result for compact surfaces X without 
boundary. If X has no boundary, let be the inner angle of the triangle AviVjVk at 
Vi is dehned to be zero if Vi,Vj,Vk do not form a triangle). Then the area Aijk of 
AviVjVk satishes the equation, 

+ 9f + = TT + 

Since ai = ^ we have 

5 ;^ K, ^2^N-'£ Of 

i i,j,k 

= 2tiN - ( 6 »f + 9f + 9f) 

AviVjVk 

= 2tiN — 7r|F| — \Area{X) 

= 2'kx{X) — XArea{X). 

As a consequence, if then the average curvature of a circle-packing metric 

is 27rx{X)/N. In particular, the constant curvature metric has curvature 27rx{X)/N. 

3.2. The Combinatorial Ricci Flow 

We dehne the combinatorial Ricci flow (with background metric AT^) to be 

(3.1) dri{t)/dt = -Kis{ri{t)). 

The evolution of curvature was hrst obtained by He [He] for zero weighted triangula¬ 
tions in Euclidean geometry. 

Proposition 3.2. Under the Ricci flow (3.1), the curvature Kiit) evolves according 
to 

dK,/dt = Y CrjiKj - Ki) + XB,K, 
j~i 
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(3.2) 



where the sum is over all vertices Vj adjacent to Vi, Cij and Bi are positive elementary 
functions in the radii ri, ...,rAr. Furthermore, Cij = Cji. 

This is a direct consequence of proposition 2.4 since dKi/dt = — dOl^/dt where 

are inner angles at Vi in the triangle AviVjVk- 

As a consequence, we obtain the maximum principle for curvature evolutions. 

Corollary 3.3. (The Maximum Principle) Let rit) = (ri(t),..., riv(t)) he a solution 
of the Ricci flow (3.1) in an interval. Let M{t) = max{iPi(t),..., iPAr(t)} and m(t) = 
min{A:i(t),..., iPjvW}- 

(1) In the case of Euclidean geometry E‘^, then the function M{t) is non-increasing 
in t and m(t) is non-decreasing in time t. 

(2) In the case of hyperbolic geometry , then max(M(t),0) is non-increasing in t 
and min(m(t), 0 ) is non-decreasing in time t. 

3.3. The long time existence of solutions 

As a consequence of the maximum principle, we have. 

Proposition 3.4. Suppose the background geometry is or . For any initial 
metric r(0) G R>o, the solution to the Ricci flow (3.1) with the given initial metric exists 
for all time t >0. 

Indeed, we note hrst of all the curvature Ki is always bounded, i.e., 

(2 — (i) 7 r < Ki < 2n 

where d is the degree of the vertex Vi. 

In the case of Euclidean geometry by (3.1), we obtain 

for some positive constants ci and C 2 . Thus rflf) remains bounded above and away from 
zero as long as t is bounded. This implies that the Ricci flow (3.1) has a solution for all 
time t > 0 for any initial value r( 0 ) G rt>o- 

In the case of since Ki(t) = d(ln{coth{ri/2)))/dt, we obtain 

cothri/2 < 

where c is some constant. Thus ri{t) remains bounded away from 0 as long as time t is 
bounded. On the other hand, the maximum principle (corollary 3.3) shows that ri{t) is 
bounded from above as well. To see this, we hrst have. 

Lemma 3.5. For any e > 0, there exists a number L so that when ri > L, the inner 
angle 6i in the hyperbolic triangle AviVjVk is smaller than e. 
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This is a simple fact from hyperbolic geometry. Indeed, if we put Vi in the Euclidean 
center of the disk model for the hyperbolic plane and choose a large radius circle centered 
at Vi in the hyperbolic space, then the Euclidean diameters of the other two circles centered 
at Vj and Vk can be made very small. Thus the angle at Vi can be made arbitrary small if 
the radius is large. 

As a consequence of the lemma, if there exists a sequence of time tn ^ a where 
a G (0, cx)] so that Tiitn) tends to inhnity, then Ki(tn) will tend to 2%. This contradicts 
the non-increasing property of the function M{t) < 27r. 

By combining the above two, we see that as long as time remains bounded, then ri{t) 
and l/ri{t) remain bounded. This implies that the solution to (3.1) exists for all time 
t > 0 . 

As a consequence of lemma 3.5 and the maximum principle, we have. 

Corollary 3.6. Suppose r(t) (t > 0) is a solution to the Ricci flow (3.1) in hyperbolic 
background geometry, then there exists a constant C so that rfli) < C for all time t > 0 
and all indices i. 

3.4. An existence theorem in the hyperbolic and Euclidean cases 

The following proposition shows that it suffices to establish a compactness result in 
order to show the long time convergence of the solution to (3.1). 

Proposition 3.7. Suppose r{t) fort G [0, cx)) is a solution to the Ricci flow (3.1) in 
geometry so that the set {r{t)\t G (0, cx)} lies in a compact region in R>o- Then r{t) 
converges exponentially fast to a circle packing metric in R>o whose curvature at each 
vertex is zero. 


Proof. Since the coefficient R in (3.2) is an elementary function in ri,..., rjv and is always 
positive, it follows from the assumption that there exist two positive constants ci and C 2 
so that 

Cl < Bflnit), ...,riv(t)) < C 2 

for all time t > 0. Let Mflf) = max(M(t),0) and m^{t) = min(m(t),0). By equation 
(3.2), we obtain 

dMflf)/dt < —ciMflf) 


and 


dmflf)/dt > —C 2 m*(t). 


Thus there are two constants C 3 > 0 and C 4 > 0 so that 


cse < m*(t) < M*(t) < c^e 


This shows that the curvature Kflf) converges exponentially fast to zero. This in turn 
implies that fg Ki(s)ds converges exponentially fast to some constant. Since (3.1) says 
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coth(ri(t)/2) = ce-Jo it follows that limt^oo= -Ri exists in R>o and ri(t) 

converges to Ri exponentially fast. Thns the resnlt follows. QED 

The same proposition also works for the Enclidean case. In this case, we shonld 
normalize the eqnation. By the Ganss-Bonnet theorem, the snm of the cnrvatnres is 
27rx(-^)- Thns the average cnrvatnre Kav is 27rxiX)/N which is independent of time t. 
We dehne the normalized Ricci flow to be 

(3.3) dri{t)/dt = -{Ki - Kav)ri 

It is easy to verify that r{t) solves (3.3) for = E‘^ if and only if e^“^*r(t) solves 

(3.3) . 

Proposition 3.8. Suppose r{t) for t G (0, cx)) is a solution to the normalized Ricci flow 

(3.3) in E‘^ so that the set {r{t)\t G (0, cx))} lies in a compact region in R>o- Then r{t) 
converges exponentially fast to a circle packing metric in R^g whose curvature at each 
vertex is the average curvature Kav 

Proof. To show that the metrics r{f) = (ri(t),..., r_v(t)) converge exponentially to a vector 
in R>o, by the same argnment as in the proof of proposition 3.7, it snfhces to show that 
the cnrvatnres K(t) = {Ki(t),KN^t)) converges exponentially fast to Rra„(l,..., 1). To 
this end, let ns consider the fnnction g{t) = —KavY- Since -R* = NKav, 

we have g{t) = ~ Thns the derivative g'{t) = 2Yf!i^iKi{t)K'flf). By 

proposition 2.4 and dKav/dt = 0, the evolntion eqnation for the cnrvatnre Ki(t) satishes 
the following, 

K'flt) = - Ka.) - {K, - Ka.)) 

j~i 

= - Ki) 

j~i 

where Cij = Cji is positive and is an elementary fnnction in ri,..., rjv- Thns 

N 

= 2 -R*)- 

i=l ir^j 

By switching the order of i, J, we also have 

N 


(3,4) = 

i~j 
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where the sum is over all pairs of adjacent vertices Vi and Vj. Since the metrics r{t) lies in 
a compact region in R >05 f^e coefficients C'ij(ri(t), ...,rN{t)) remain bounded for all time 
t > 0, i.e., Cij > Cl > 0 for some constant ci. 

Claim. There exists a constant C 2 > 0 depending only on the triangulation T so that 

N 

C 2 < Y, 

i=l i~j 

for all time t > 0. 

Assuming the claim, by (3.4), we obtain the following 

g'{t) < -C 2 g{t). 

This implies that g{t) < for all time t > 0 for some constant C 3 . In particular, we 

obtain \Ki{t) — Kav\‘^ < gif) < Thus the proposition holds. 

To see the claim, we hrst note that since Cij are uniformly bounded from below, it 
suffices to prove that 

N 

(3.5) 5^(K, - < C4 5^(K* - 

i=l i^j 

for some constant C 4 . By Cauchy’s inequality that for any sequence {ui} 

a- ^ 

fi - 

i=i 

we obtain 

N N 

YiKr - Kav)^ < l/N Y 

i=l i,j=l 

On the other hand, since the surface X is connected, for any two vertices Vi and Vj in T, 
there exists a sequence of vertices Vm^ — Vi,Vm 2 j ■■■Xmi — Vj so that Vm,, and Vm^+i are 
adjacent. We can estimate the value (Ki — Kj)"^ in terms of the sum of adjacent terms 
{Kruh ~ ■ Tffis shows that there exists a constant C 5 depending only on T so that 

(iC, - K,)2 < C 5 - Km)^- ThuS 

N 

Y {K, - Kj)^ < c,NY{K^ - 

i,j=l i^j 

Combining all these, we establish the claim. Thus the proposition follows. 

3.5. The Ricci flow is variational 
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In this section, we establish the fact that the combinatorial Ricci flow is essentially 
the gradient flow of a fnnction. In the case that the backgronnd geometry is or H^, 
the f un ction is convex. 

Let Ui be a fnnction of ri so that dui/dri = l/s{ri). To be more precise, we choose 
Ui = Inri when Ui = lntanh(ri/2) in and Ui = lntan(ri/2) when 

Let U be the following open set. If then U = R^; if then 

U = (—00,0)^. In the case of the spherical geometry, we take U to be the inverse image 
of all possible r = (ri,..., r_v) nnder the map {ui ,..., un) ^ (ri,..., tn)- 

Under this change of variable, the Ricci flow (3.1) takes the following form 


(3.6) dui/dt = —Ki. 

Indeed, (3.6) is a gradient flow in the variable u. To see this, by le mm a 2.3, we 
have dKi/drjs{rj) = dKj/dris{ri). This is eqnivalent to dKi/duj = dKj/dui. Thns 
the 1-form Kidui is closed in the simply connected space U. Dehne a fnnction 

f{u) = Kidui where a is any point in U. Then we have df{u)/dui = Ki. This 

shows that (3.6) is the negative gradient flow of the f un ction /. We remark that a family 
of fnnctionals like this / was hrst constrncted by Cohn de Verdiere in his paper [CV] in 
the case the weight fnnction is zero. Cohn de Verdiere also hrst proved the convexity of / 
in the zero weight case. 

Proposition 3.9. (a) If , then the function f{u) : U 'R. is strictly convex. 

(b) If then the function f{u) : U —R is convex. Furthermore, it satisfies 

f(u -|- a(l, 1 ,..., 1 )) = f{u) for all a and f is strictly convex when restricted in the plane 

P = {(wi, ...,'Ujv)| = 0}- 

Proof. Let ns consider the hessian of the f un ction /. 

Let aij = d'^f/duiduj = dKi/duj. By lemma 2.2, we have an > 0 and aij < 0 ( 
Oij < 0 if is adjacent to Vj). Fnrthermore, lemmas 2.2 and 2.3 show that if 
then d6i/drismh.{ri) +d6i/drjsmh{rj) + d6i/drkimah.{rk) < 0. Now snm over all snch 
ineqnalities at the i-th vertex and nse the fact that a^- = —^^^d6l^/drjsmh.{rj), we 

obtain > 0- Recall a simple fact from linear algebra. 

Lemma 3.10. Suppose A = [aij]nxn is a symmetric metrix. 

\aij\ for all indices i, then A is positive definite 
(b) If an > 0 and aij < 0 for all i ^ j so that — 0 U then A is 

semi-positive definite so that its kernel is one-dimensional. 

Thns part (a) of the proposition follows. 

To see part (b) of the proposition, note that lemmas 2.2 and 2.3 show that an > 0, 
aij < 0 and Ei=i ~ 0 - linear algebra fact above, this implies that the matrix 

[aij] is semi-dehnite. Thns / is convex in R"". The identity /('u-|-(a, a,..., a)) = f{u) follows 
from the fact that Ki{ar) = Ki{r). The last fact that /|p is strictly convex follows from 
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the linear algebra consideration, namely, the kernel of the matrix [aij]NxN is generated by 
the vector (l,l,...l). 

To prove the lemma 3.10 in the case of (a), it snfhces to show that det(A) > 0, since 
the same argnment will show that all determinants of the principle snbmatrices of A are 
positive. This in tnrn implies that A is positive dehnite. Now to show det(A) > 0, we 
hrst show that det(A) ^ 0. This follows easily from the given condition. For instance, 
if {xi,...,Xn) were a non-zero vector so that Ylj^ijXj = 0 for all i, take Xk to be the 
coordinate with the maximnm absolnte valne. Now consider O'kjXj = 0. We obtain 

ttkk\xk\ = I '^ij^k^kiXil- Bnt this contradicts the assnmptions that akk > \^ki\ and 

> nraxillxi|}. Next mnltiply all off diagonal entries of A by a constant t so that 
\t\ < 1 to get a new matrix B{t). The matrix B{t) satishes the same conditions that 
A satishes, thns det(i?(t)) 7 ^ 0. Bnt det(i?(0)) > 0, thns det(i?(t)) > 0. In particnlar 
det(A) = det(i?(l)) > 0 . 

To see part (b) of the lemma, hrst of all it is clear that [aij] is semi-positive dehnite 
(since if we replace an by an -t-s for a small positive nnmber s, it becomes positive dehnite 
by part (a)). Thns all eigenvalnes of it are non-negative. Now it is easy to see that the 
kernel is one-dimensional spanned by (1,1,, 1). Thns the resnlt follows. QED 

Given a smooth strictly convex fnnction dehned on an open convex set H, the gradient 
map sending a point in D to the gradient of / at the point is always injective. Thns as a 
conseqnence of the above proposition, we obtained a new proof of the nniqneness part of 
Thnrston’s theorem. 

Corollary 3.11. ([Th]) Let n : R>o be the map sending a metric (ri,...,rjv) to 

the corresponding curvature (iFi, ...,Kn) where or E‘^. Then 

(a) in the hyperbolic background geometry, n is injective, i.e., the metric is determined 
by its curvature. 

(b) in the Euclidean geometry, n restricted to the subset {(ri,..., rjv) G R>o|n^]^ri = 
1 } is injective, i.e., the metric is determined by its curvature up to a scalar multiplication. 

§4. Degeneration of Circle Packing Metrics 

We will recall a related work of [Th]. For completeness a slightly diherent proof is 
presented. The proof below follows closely [MR]. For simplicity, if / = {r^, ...ri} is a snbset 
of vertices V, we will also nse I to denote the index set {j, 

Proposition. 4.1 ([Th]) Suppose (T, $) is a weighted generalized triangulation of a closed 
surface X and I is a subset of vertices V. Here the weight is a map $ : i? —[ 0 , 7 r). Let 
rl"") = {r^\ ..., be a sequence of circle packing metrics based on (T, $) in background 
geometry so that lim^^oo ~ 0 for i E I and lim^^oo rl^'^ > 0 for i ^ I. Then 

Urn ^Ki{r^^^) = - (tt - $(e))+27rx(F/), 

i^I {e,v)^Lk{I) 
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where Fj is the CW-subcomplex consisting of cells whose vertices are in I and Lk{I) = 
{(e, v)\e is an edge so that e fl / = 0 and the vertex v E I so that e, v form a triangle}. 

Furthermore, if the weight $ : i? —> [0, 7 r/ 2 ] and the background geometry is or , 
then for any circle packing metric r based on (T, $) and any proper subset I of vertices V, 
we have 

'^Ki{r)>- (7r-$(e)) + 27rx(i^/). 

iel {e,v)eLk{I) 


Proof. Recall that 9{^ denotes the inner angle at Vi in the triangle AviVjVk- Let X ^ X 
be the nniversal cover. Consider triangles in T having a vertex in /. These triangles can 
be classihed into three types Ai, A 2 and A 3 where a triangle A is in Ai if and only if a 
lift of it in the nniversal cover has exactly i many vertices in 7 r“^(/). (Note that in the 
snrface X, there may be triangles having all vertices being Vi.) 

Now the cnrvatnre Ki = 2 ti — where ai is the cone angle at Vi. Here Yhiei 
be written as 

E E (9f+»f)+ E (»/*+»?+«?)■ 

Vi^I ^AviVjVk^Ax Vi,Vj^I,AviVjVk^A2 AviVjVk^As 

(n) 

It is known that as the radins r} ^ tends to zero, the following holds. 

(1) If AviVjVk E A 3 , then lim^^oo(^i^(?'*'"'^) + = n. Indeed, the 

triangle is shrinking to a point in and thns is becoming a Enclidean triangle. 

(2) If AviVjVk E A 2 and i,j E I, then = tt. Indeed, the 

triangle degenerates to a geodesic segment. 

(3) If AviVjVk E Ai and i E I, then lim^^cx) = tt — ^{vjVk) by dehnition. 

Here VjVk E Lk{I). 

To snmmary, we obtain 


lim = - Y. (tt - $(e)) + 27r(|/| - I/ 2 IH 2 I - I/ 2 IH 3 I). 

iel {e,v)eLk{I) 

Here \Ai\ denotes the nnmber of triangles of type Ai. Note IH 3 I is the nnmber of 
2-cells in Fi and |/| is the nnmber of 0-cells in Fj. By the constrnction, we see that the 
nnmber of edges in Fj is 1 / 2 (|H 2 | + 3 |H 3 |). Thns |/| — I/ 2 IH 2 I — I/ 2 IH 3 I is the Enler 
characteristic of Fj. This ends the proof. 

To see the second part of the proposition, we note that since the backgronnd geometry 
is or the snm of the inner angles of a triangle is at most tt. Thns, in the hrst 
limit (1) above, tt is the maximal valne. By the same argnment, in the second limit 
( 2 ) above, tt is again the maximal valne which is never achieved by a non-degenerate 
triangle. In the third limit (3), since $(e) < 7 r/ 2 , nsing lemma 2.1 on the monotonicity, 
we see that tt — ^{vjVk) is also the maximal valne for all possible circle packing metrics 
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based on (T, $). Furthermore, since I ^ V, the cases (2) or (3) do occur. This shows 
that —Yli{e v)eLk{i)^'^ ~ ^(^)) + 27rx(Fj) is strictly less than the sum of the curvature 
Ki{r) for any circle packing metric based on (T, $). QED. 

Remark. The proof works for all CW -cell decompositions so that each cell has 3 vertices 
and each vertex has degree at least 2 . 

Corollary 4.2. Make the same assumption as in proposition 4-1 assume that the 
weight ^ : E ^ [0,7r/2]. Suppose - E(e,^;)€Lfc(/)(^ ~ ^(^)) +^^x{Fi) > 0. Then, 

(a) either there exists a null homotopic loop consisting of 3 edges 61 , 62,63 so that 
Ell $( Ci) > TT. Furthermore, if X ^ 3“^, then these 3 edges do not form the boundary of 
a triangle, or 

(b) X = S'^ and I — V — {n} for one vertex v. 

Indeed, choose I' be the non-empty subset of / of the smallest cardinality so that 
~ E(e v)eLk{i)('^ ~ ^( 6 )) + 27rx(F7) > 0 . In this case, F// is connected since the quantity 
is additive with respect to connected components of Fj. Thus x(-^/0 < 1 so that the 
equality holds if and only if Fj/ is contractible. Since v)eLk{i')i'^ ~ ^( 6 )) — O 5 

see that xi.T'r) = 1 and Fj/ is contractible. Let 61 ,..., 6 fc be the set of all edges in Lk{F). 
They form a null homotopic loop due to the contractibility of Fj/. Now if Lk{I') 7 ^ 0, 
then v)eLk{i')('^ ~ ^(6)) — kj27r. Thus k = 3. Note k = 1,2 are ruled out by the 
assumption on the generalized triangulation T. If these three edges form the boundary of 
a triangle, then Fp must be the complement of the triangle. Thus X = S'^. 

If Lk{I') = 0, then V — F must consists of vertices so that no two of these vertices are 
adjacent. On the other hand, F// is contractible. Thus V — F consists of only one point. 
Thus we see that X is a union of two contractible sets and X = S'^. 

Corollary 4.3. Make the same assumption as in proposition 4-1 assume that the 
weight $ : F —> [0, 7 r/ 2 ]. If there are no subsets of V satisfying the condition in corollary 
4.2 and if -J2{e,v)eLk{i)i^ “ ^(e)) + 27rx(F/) = 0, then 

(a) either there exists a null homotopic loop consisting of k = 3 or 4 edges 61 , ...,ek 

so that Ef=i = {k — 2 )ti. Furthermore, if X ^ S"^, then these k edges do not form 

the boundary of a triangle, or union of two adjacent triangles, or 

(b) X = S‘^ and I = V — {v,v’\ for two non-adjacent vertices v,v', or 

(c) X = RF2 andI = V - {n}. 

Let I' be the subset of the smallest cardinality so that E(e v) eLk{i)i^ - ^(e)) + 
2tix{Fi) = 0. By the assumption, we conclude that F// is connected. Thus x{Fi/) < 1. 
By the hypothesis, x{Fii) > 0. If x{Fr) = 1, then by the given hypothesis, Fj/ 7 ^ 0. By 
the same argument as above, we obtain conclusion (a). If xiFj) = 0, then Lk{F) = 0, 
i.e., V — F consists of vertices so that no two of them are adjacent. But Fj/ is either the 
annulus or the Mobuis band. This shows that either X = S‘^ and F = V — {v,v'} where 
V, v' are not adjacent or X = RF^ and F ~V — {n}. 
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Corollary 4.4. Suppose X is a closed surface of negative Euler characteristic and 
(T, $) is a generalized triangulation on X so that the conditions (1.5a) and (1.5b) in theo¬ 
rem 1.2 hold. If r^'^^ G R>q is a sequence of circle packing metrics in the hyperbolic back¬ 
ground metric based on (T, $) so that > 0 and limsup^^j^ Kifr^'^^) < 

2ti for all indices i, then contains a convergent subsequence in R>o- 

Indeed, by lemma 3.5 and limsup^^o^ < 2ti, we see that r^f^^ is bounded 

from above. To see that is bounded away from 0, suppose otherwise that it contains 

a subsequence converging to zero. Let I be the non-empty subset of indices so that 

limj^oo= 0 for t G / and lim^^oo> 0 for f ^ /. Since lim inf > 0 

for all indices i, we obtain lim^^oo 'Yhiei > 0. By corollaries 4.2 and 4.3, we see 

that one of the conditions (1.5a), (1.5b) in theorem 1.2 must be false. This contradicts the 
assumption. 

§5. A Proof of Theorem 1.2 

We will prove the following result which is slightly more general than theorem 1.2. 

Theorem 5.1. For any weighted generalized triangulation (T, $) of a closed connected 
surface X of negative Euler characteristic and any initial choice of circle-packing metric, 
the solution to (l.f) with the given initial value exists for all time and converges if and 
only if the following two conditions hold. 

(1.5a) for any three edges ei, 62 , 63 forming a null homotopic loop in X, if'Ylii=i ^(^i) > tt, 
then these three edges form the boundary of a triangle of T, and 

(1.5b) for any four edges ei, 62 , 63 , 64 forming a null homotopic loop in X, — 

2 %, then these Ci’s form the boundary of the union of two adjacent triangles. 

Furthermore, if it converges, then it converges exponentially fast to a hyperbolic metric 
on X without any singularity. 

Note that if the weight $ is zero, then conditions (1.5a) and (1.5b) hold automati¬ 
cally. Thus the above result implies the result of Andreev-Koebe-Thurston that for any 
triangulation of a closed surface of negative Euler characteristic, there exists a hyperbolic 
metric on the surface and a tangential circle packing in the metric so that the combinatorial 
pattern of tangency is given by the triangulation. 

The proof of the theorem goes in two steps. In the hrst step, we give a new proof of 
Thurston’s existence of circle packing metric using Ricci flow. This implies that the strictly 
convex function / whose negative gradient is dui/dt = —Ki where Ui = lntanh(ri/2) has 
a critical point. In particular, this implies that the function / is bounded from below. On 
the other hand, corollary 4.4 shows that the function / satishes the Palais-Smale condition. 
Thus the negative gradient flow of / converges to the critical point of /. Combining this 
with proposition 3.6, we see that all gradient lines converge exponentially fast to the critical 
point. 


17 



5.1. In this subsection, we produce a new proof of the existence of the zero curvature 
circle packing metric based on the weighted generalized triangulation (T, $). 

To begin, let us choose the initial metric r(0) = (ri(0), ...,r7v(0)) with all ri(0) suffi¬ 
ciently large so that it curvature at Vi is close to 27r. This is possible by lemma 3.5. In 
particular, we have /ci(0) > 0. Let r{t) = (ri(t),..., rNif)) (t > 0) be a solution to the Ricci 
flow (3.1) in background metric with the initial value r(0). Then by the maximum 
principle corollary 3.3, we have ki{t) > 0 for all time t > 0. On the other hand, by the 
maximum principle again we have limsup^^oo/ci(t) < 27r. Thus by corollary 4.4, we see 
that the set {r(t)\t G [0, cx))} lies in a compact region in R>o- Thus by proposition 3.7, 
the solution r{t) converges exponentially fast to the zero curvature metric. In particular, 
this shows the zero curvature metric exists. 

5.2. To show the convergence of the Ricci flow (3.1) in background metric in the 
general case, we use the result in subsection 3.5. Recall that if we consider the change of 
variable Ui = lntanh(ri/2), then the Ricci flow becomes dui/dt = —Ki which is proved to 
be a negative gradient flow. Let / be a potential function. It is known by proposition 3.9 
that the function / is strictly convex dehned on (— cxd, 0)^. We claim that / satishes the 
Palais-Smale condition. As a consequence, we see all gradient lines converge to the critical 
point. Thus the result follows. 

To verify the Palais-Smale condition, hrst of all by §5.1, we see that / has a critical 
point which must be the minimal point due to the convexity. Thus / is bounded from 
below. To verify the second condition in Palais-Smale, take a sequence of points } in 
(—CX), 0)-^ so that the gradient of / at converges to zero. Let the corresponding points 
of be G R>o- Then we have that the curvature tends to zero for all 

indices i as n tends to inhnity. By corollary 4.4, this implies that the set e Z>o} 

lies in a compact set. This in turn implies that contains a convergent subsequence 

in (—cx,0)'^. This ends the proof. 

Finally the necessity of the conditions (1.5a) and (1.5b) which was proved in [Th] 
follows from proposition 4.1. Indeed, since the zero curvature metric exists, we conclude 
that for any subset / G R, the quantity — X](e i;) eLk{i)i^ - ^(e)) + 27 rx(Fj) < 0. On the 
other hand, if (1.5a) or (1.5b) is false, then take / to be the vertices in the disk bounded 
by the null homotopic loop formed by the edges. We see that —J2{ev) eLk{i)i^ - + 

27r(x(F/)) > 0. 

Remark 5.3. The proof in fact shows the following for a generalization of the Ricci flow 
in hyperbolic background metric. Namely, suppose there exists a circle packing metric R 
so that it curvature is Ci at the i-th vertex. Then the following flow dri/dt = —{Ki — 
Ci) sinh(ri) has a solution for all time t > 0 for all initial metrics and the solution converges 
exponentially fast to R. 

§5. A Proof of Theorem 1.1 and Proposition 1.3 

We will prove the following theorem for generalized triangulations which implies the¬ 
orem 1.1. 
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Theorem 6.1. For any weighted generalized triangulation T of a closed connected surface 
X and for any initial circle-packing metric based on the weighted triangulation, the solution 
to the normalized combinatorial Ricci flow (1.2) in Euclidean background metric with the 
given initial value converges if and only if for any proper subset I C V , 


(1.3) 27rx(X)|/|/iV>- (7r-$(e)) + 27rx(Fz). 

{e,v)ELk{I) 

Furthermore, if the solution converges, it converges exponentially fast to the metric of 
constant curvature 2 txx{X)/N. 

The proof is essentially the same as in §5 and goes in two steps. In the first step, we 
show that the constant curvatnre circle packing metric exists dne to condition (1.3). Next 
we show that the normalized Ricci flow is the negative gradient flow of a convex fnnction 
in the variable Ui = Inr^ and that the fnnction / satisfies the Palais-Smale condition. 
Combining this with proposition 3.6, we see that all gradient lines converge exponentially 
fast to the critical point. 

6.1. In this snbsection, we prodnce a proof of the existence of the constant cnrvatnre 

circle packing metric based on the weighted generalized triangnlation (T, $) nnder condi¬ 
tion (1.3). This resnlt is implicitly in [Th]. Onr proof follows the work [MR]. Let P be 
the set {r = (ri, ...,r_;v) G R>ol = !}• Consider the cnrvatnre map 11 : P —> 

sending the metric r to its cnrvatnre n(r) = {Ki{r),KN{r)). By corollary 3.9, it is 
injective. Fnrthermore, by Ganss-Bonnet, its image lies in the hyperplane defined by the 
eqnation — 27rx(-T). Consider the convex polytope Y in the hyperplane defined 

by the set of all linear ineqnalities y Xi P 'Tjiei (tt — $(e)) -I- 2tix{Fi) where / is a 
proper snbset of vertices. Then by proposition 4.1 the map II : P sends compact sets 
to compact sets. Since both P and Y are homeomorphic to R^“^, by the invariance of 
domain theorem, the map 11 : P —T is a homeomorphism. In particnlar, we see that the 
constant cnrvatnre metric exists if and only if (1.3) holds for all proper snbsets / of V. 

6.2. To show the convergence of the Ricci flow for all initial choice of metrics, let ns 

take Ui = Inr^. Then the normalized Ricci flow becomes dui/dt = —{Ki — Kav)- This later 
eqnation is the negative gradient flow of a convex fnnction f{u) defined on the n-space R^. 
Fnrthermore, if we restrict to the snbspace Z = {u = (wi, ...,un) G R^^I = 0}? f^e 

fnnction / is strictly convex proposition 3.8. By the resnlt in snbsection 6.1, we see that / 
has a critical point. Thns it is bonnded from below. We claim that / satisfies the Palais- 
Smale condition on the plane Z. To this end, take a seqnence of vectors in Z so that 
the gradient of / at tends to zero. Let r^^'> be the corresponding points in the space 
P. Then lim^^oo = Kav for all indices i. By the assnmption that Kav{l, ■■■, 1) lies 

in the open cell n(P), we see that lies in a compact set in P. Thns contains a 

convergent snbseqnence. This shows all gradient lines converge to the constant cnrvatnre 
metric. Thns by proposition 3.6, the solntion r{f) converges exponentially fast to the 
constant cnrvatnre metric. 
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6.3. The proof in fact shows the following for a generalization of the Ricci flow in 
Enclidean backgronnd metric. Namely, snppose there exists a circle packing metric R so 
that it cnrvatnre is Ci at the i-th vertex. Then the following flow dri/dt = —{Ki — Ci)ri has 
a solntion for all time t > 0 for all initial metrics and the solntion converges exponentially 
fast to R. 

6.4. A proof of proposition 1.3 

Proposition 1.3 (a) If the Euler characteristic of the surface X is positive, then condition 
(1.3) is valid for all generalized triangulations of X so that there are no null homotopic 
loops consisting of 3 edges ei,..., 63 so that Ylii=i ^ 

(b) If the Euler characteristic of the surface X is zero, then (1.3) is valid for all 
triangulations of X so that (1.5a) and (1.5b) hold. 

This follows from corollaries 4.2 and 4.3. 

§7. Some Questions 

There are several questions which we hnd interesting relating to the results in the 
paper. 

1. Generalizing Lemma 2.2 

It can be shown that lemma 2.2 is false if the assumption G [0,7r/2] is relaxed to 
$i’s in [0, tt). In fact there exists a choice of G [0? ^r) so that ddijdri < 0. But on 
the other hand, it seems for applications, the convexity of the function / in proposition 
3.9 is essential. Thus the question is for what choices of $n’s in [0, tt) is the matrix 
[s{rj)d6i/drj\^xz positive or semi-positive dehnite. 

2. 2-Dimensional Combinatorial Ricci flow for all piecewise constant curvature met¬ 
rics 

Hamilton’s ricci flow works for all Riemannian metrics. Fix a triangulation on a 
surface, the analog of Riemannian metric seems to be the piecewise constant curvature 
metric obtained by assigning lengths to each edge in the triangulation. Is there any Ricci 
flow dehned on the space of all these metrics? This may have some applications to the 
problems like what is the space of all geometric triangulations of a triangle with the hxed 
combinatorial data. 

3. Investigate the 2-dimensional combinatorial Ricci flow in theorems 1.1 and 1.2 so 
that the combinatorial conditions (1.3) and (1.5a) and (1.6b) are not valid. The solution 
still exists for all time t > 0. Is it true that for any solution r(t) = (ri(t),..., rAr(t)) where 
t >0, the limit lim^^oo fi{t) always exists in the extended set [0, cx)]? 

In the case of the Euclidean background metric, for negative Euler characteristic closed 
surfaces, if r{t) with t > 0 is a solution to the Ricci flow (1.3), is it true that the metric 
r{t) will have non-positive curvature when t is large? 
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4. Investigate the 2-dimensional combinatorial Ricci flow in the case that the weight 
is in [ 0 , 7 r). For instance, can one prodnce a new proof of the resnlts of Igor Rivin [Ri] 
nsing heat eqnations? 

5. We showed that rigidity of the circle packing follows from the strictly convexity of 
the potential fnnction in this paper (corollary 3.11). The convexity is based on the main 
technical lemma 2.3. Derive the main resnlt of He [He] on the rigidity of locally finite 
infinite circle packings from the convexity of the potential fnnction. 

6 . Is there any 3-dimensional combinatorial analogons of Hamilton’s 3-dimensional 
Ricci flow [Ha2]? 

The recent work of David Glickenstein’s UCSD thesis work ([Gl]) seems to indicate 
that the analog of the Yamabe flow exists in the combinatorial setting in dimension 3 (see 
also [GR]). 

Appendix. Some computations involving triangles 

Al. Basic results on triangles 

Given a triangle AviVjVk in one of the three geometries and H^, let 

the inner angle at the Vi vertex be 6 i and the length of the edge VjVk be Xi. Recall that 
the fnnction s(a:) is defined according to the geometry as follows: s{x) — x for E‘^, 
s{x) = sin(a:) for 3“^ and s{x) = sinh(a:) for Let Aijk — s{xi)s{xj) sm{9k)- By the 
sine law, the qnantity Aijk is symmetric in i,j, k. The following is a very nsefnl lemma for 
compnting derivatives. 

Lemma A-1. In all three geometries we have, 

(a) dOijdxi — s{xi)/Aijk- 

(b) dOi/dxj = —d9i/dxiCos{9k)- 

(c) d9ildxjs{xi) is symmetric in i,j. 

(d) Consider Xk as a function of Xi,Xj,9i. Then 

dxkjdxi = cos(6*j). 


The proof is a simple application of the cosine law. We will omit the details. It is 
interesting to note that these derivatives formnlas are the same for all three geometries. 
One interesting conseqnence of lemma A-l(c) is the following. 

Corollary A-2. For any triangle in the geometry the 1-form j k} ^i^{^i)dxi in 

variables Xi,Xj,Xk is closed. 

A 2 . A proof of lemma 1.3 

A simpler proof nsing mathematical software Maple VH is given in the next snbsection. 
In this snbsection, we present a conventional proof. 
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Suppose we have a collection of three closed disks Ci, Cj, and Ck centered at Vi, Vj and 
Vk so that their radii are r^, rj and rk and their exterior intersection angles are (pi, (pj, pk ^ 
[0, tt) in geometry K^. Consider the triangle formed by Vi,Vj,Vk- Let its inner angle at 
Vn be 9n and the edge length opposite to Vn be Xn- We want to show that d6i/drjs{rj) is 
symmetric in i,j. 

First of all if Li, Lj and Lk are the geodesic lines passing through the pairs of the 
intersection points of {Ck,Cj}, {Ck,Ci} and {Ci,Cj} {Li is the tangent to Cj and Ck if 
they intersect in one point), then the three lines Li, Lj, Lk intersect in one point. 

To see this, we note that the result is well known in the Euclidean geometry To 
see the result in hyperbolic geometry, let O be the point of intersection of Li and Lj. 
Take the disk model for so that O is the Euclidean center. Then both Lj and Lj are 
Euclidean line segments. Let L). be the Euclidean diameter passing through O and one 
of the intersection point of Ci with Cj. The by the result in the Euclidean case, we see 
that Lk must contain the other point of intersection of Ci and Cj. This shows, L'j^ = Lk. 
The same argument also works for the spherical case where we consider the stereographic 
projection from the antipodal point of O. Let Uj. be an intersection point of the circles Cg 
and Ct where {r, s, t} — {i,j, k}. 

To begin the proof, by lemma Al, we have 

s{rj)d9ijdrj = s{r j){d9ij dxidxij dr j + dOijdxkdxkjdrj) 

= s{rj)d6i{dxi{dxijdrj — cos{6 j)dxk / dr j) 

^ ^(xj) sm{9j) s{rj){cos{ef) - cosjOj) cosjef)) 

Aijk sin{9j) 

Here 6**^ and 9^ are the inner angles at Vj in the triangles AvjViUk and AvjVkUi. 
Note that the hrst fraction s{xi) sin(^j )/Aijk is symmetric in i, j due to the sine law. Thus 
it suffices to show that 

s{rj){cos{9y) — cos{9j) cos{9j^')) 
sin{9j) 

is symmetric in i,j. 

To this end, we will need the following lemma about quadrilaterals in so that two 
opposite angles are right angles. 

Lemma A2. Suppose Q is an immersed quadrilateral in geometry K‘^ so that its four 
edges are T,X,Y,Z counted clockwise. Suppose the inner angle between sides {T,X} and 
{Y, Z} are right angles and the inner angle between sides {X, Y} is 9, and edges T, Y are 
disjoint. Let the lengths of sides X, Y, Z be x, y, z and t be the length of T of X f] Z = ^ 
and t be the negative of the length of T if X intersects Z. Then the following holds. 

(a). In the Euclidean geometry E‘^, 

t = {y — X cos 9) / sin 9. 
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(b) . In the hyperbolic geometry , 

tanh(t) = cosh(a:)(tanh(y) — tanh(a:) cos( 6 *))/ sin( 6 *). 

(c) . In the spherical geometry S'^, 

tan(t) = cos(a:)(tan(y) — tan(a:) cos( 6 *))/ sin( 6 *). 


The proof is a direct application of trigonometries for right angled triangles and qnadri- 
laterals with three right angles. (See for instance [Vi]). 

Another simple fact we shall nse is the following formnlas for right angled triangles. 
Snppose a right angled triangle has three edge lengths a, 6 , c and whose opposite angles 
are a,/3 and tx/2. Then we have, cos(a) = b/c in cos(a) = tanh( 6 )/tanh(c) in 
and cos(a) = tan( 6 )/ tan(c) in 

Let the distance from the common intersection O to the intersection point of the 
lines Lfc and ViVj be d. Let T> be d if O is inside the triangle and —d otherwise. It is 
symmetric in i,j. Then by the above le mm a and the trigonometric formnlas, we obtain 
that the expression s(rj)(cos(dj^* ) — cos(dj) cos(d*^ ))/siniOj) is eqnal to the T>, tanh(T>) 
and tan(T>) respectively. This verihes the symmetry. Note that if all angles are in 
[ 0 , 7 r/ 2 ], then the common intersection point O is inside the triangle. Thns d > 0 which is 
what lemma 2 . 2 (b) states. 

A3. A proof of lemma 2.3 using Maple VII 

The simplest proof of the lemma is to nse some mathematical software so that we do 
not have to go over all possible conhgnrations. We will present snch a proof nsing Maple 
VII below. 

Let T 3 be the space of all possible conhgnrations of pairwise intersecting three circles in 
np to isometry so that their centers form a triangle. The we claim that T 3 is connected. 
To see this, since the space of all triangles np to isometry is connected, it snffices to show 
that the snbspace of all pairwise intersecting three circles with hxed centers is connected. 
To see the later space, let ns exam the three radii. We can choose the radins of the hrst 
circle arbitrary as long as it lies in its natnral interval (e.g, it is (0, cxd) for and and 
is ( 0 , 7 r/ 2 ) for S‘^). Now the radins of the second circle mnst also be in an interval since 
it intersects the hrst circle. The radins of the third mnst be in an interval as well since it 
intersects the hrst and the second circle and these two circles intersect. This shows that 
the space is connected. 

Now to show d9i/drjs{rj) is symmetric in i,j, it is mnch convenient for the machine 
to verify that {d6i/drjs{rj)Y is symmetric in f,j. On the other hand, both fnnctions 
d6i/drjs{rj) and ddjjdrisiri) are analytic in variables $fc,ri,rj,rfc and dehned on 

a connected space so that they coincide in the snbspace = r^. Thns d6i/drjs{rj) mnst 
be symmetric in i, j if its sqnare is symmetric. 

To verify the symmetry of {{d9i/drj)rj)‘^ in i,j for Enclidean geometry E‘^ nsing 
Maple VII, we rnn the following. 
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> yi := sqrtirp + rk'^ + 2 * rj * rk * bi); 

> yj := sqrt{rP + rk‘^ + 2 * ri * rk * bj); 

> yk := sqrtirp + ri^ + 2* rj * ri * bk)-, 

> ai := arccos{{{yjf + {ykf - {yif)/{2 * {yj) * {yk)))] 

> fij ■= diff{ai,rj)*rj] 

> gl := simplify{fij)] 

> ^2 := gl^] 

> g3 := subs{{ri — rj, rj = ri, bi = bj, bj = bi}, g2) — g2-, 

> simplify{g3)-, 0 

To verify the symmetry of {{dOi/drj) smh{rj))‘^ in i,j for hyperbolic geometry 
using Maple VII, we run the following, 

> yi := arccosh{cosh{rj) * cosh{rk) + sinh{rj) * sinh{rk) * bi); 

> yj := arccosh{cosh{rk) * cosh{ri) + sinh{ri) * sinh{rk) *bj); 

> yk := arccosh{cosh{rj) * cosh{ri) + sinh{rj) * sinh{ri) * bk); 

> ai := arccos{{cosh{yj) * cosh{yk) — cosh{yi))/{sinh{yj) * sinh{yk))); 

> fij ■■= diff{ai, rj) * sinh{rj); 

> gl := simplify{fij); 

> ^2 := gl^; 

> g3 := subs{{ri = rj, rj = ri, bi = bj, bj = hi}, g2) — g2; 

> simplify{g3); 0 

To verify the symmetry of {{d9i/drj)rj)‘^ in i,j for spherical geometry S'^ using Maple 
VII, we run the following, 

> yi := arccos{cos{rj) * cos{rk) + sin{rj) * sin{rk) * bi); 

> yj := arccos{cos{rk) * cos{ri) + sin{ri) * sin{rk) * bj); 

> yk := arccos{cos{rj) * cos{ri) + sin{rj) * sin{ri) * bk); 

> ai := arccos{{—cos{yj) * cos{yk) + cos{yi))/{sin{yj) * sin{yk))); 

> fij ■■= diff{ai, rj) * sin{rj); 

> gl := simplify{fij); 

> ^2 := gl"^; 

> g3 := subs{{ri — rj, rj = ri, bi = bj, bj = bi}, g2) — g2; 

> simplify{g3); 0 

This ends the proof. 
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